Spontaneous rotating vortex lattices in a pumped decaying condensate 
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Injection and decay of particles in an inhomogeneous quantum condensate can significantly change 
its behaviour. We model trapped, pumped, decaying condensates by a complex Gross-Pitaevskii 
equation and analyse the density and currents in the steady state. With homogeneous pumping, 
rotationally symmetric solutions are unstable. Stability may be restored by a finite pumping spot. 
However if the pumping spot is larger than the Thomas-Fermi cloud radius, then rotationally sym- 
metric solutions are replaced by solutions with spontaneous arrays of vortices. These vortex arrays 
arise without any rotation of the trap, spontaneously breaking rotational symmetry. 
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While much of the possible physics of quantum con- 
densates has been examined in experiments on atomic 
gases, superfluid Helium and superconductors, there has 
recently been much interest in examples of condensates of 
quasiparticle excitations, such as excitons N, 2| (bound 
electron- hole pairs), exciton-polaritons [1, 0, 5| (super- 
positions of quantum well excitons and microcavity pho- 
tons) , and magnons (spin- wave excitations) both in mag- 
netic insulating crystals |6|, 17[ [33[ and in superfluid 
^He [1, [l, [lo|. One particular difference shown by these 
systems is that the quasiparticles have finite lifetimes, 
and as a result, they can be made to form condensates 
out of equilibrium, which are best understood as a steady 
state balance between pumping and decay, rather than 
true thermal equilibrium. The effects of pumping and 
decay in these condensates have been the subject of sev- 
eral recent works [i, [HI, [H, [H, Q E, H EM H 
which have shown that even when collisions can rapidly 
thermalise the energy distribution of a system, there may 
yet be noticeable effects associated with the energy scale 
introduced by the pumping and decay. 

The Gross-Pitaevskii equation (GPE) has been applied 
to successfully describe many features of equilibrium con- 
densates when far in the condensed regime, including 
density profiles, the dynamics of vortices, hydrodynamic 
modes — see e.g. [21] and Refs. therein. Using a mean- 
field description of the condensate, e.g. 0,0,0, one 
can recover a complex Gross-Pitaevskii equation (cGPE), 
including terms representing gain, loss and an external 
trapping potential. This letter studies the interplay be- 
tween pumping and decay and the external trapping po- 
tential in the context of the cGPE in order to illustrate 
some of the differences between equilibrium and non- 
equilibrium condensates. In the absence of trapping, 
this is the celebrated complex Ginzburg-Landau equa- 
tion that describes a vast variety of phenomena [22] from 
nonlinear waves to second-order phase transitions, from 
superconductivity to liquid crystals and cosmic strings 
and binary fluids [23] . What is of interest in this letter is 
how pumping and decay, described in the cGPE modify 



behaviour compared to the regular GPE as is widely ap- 
plied to spatially inhomogeneous equilibrium quantum 
condensates [21]. Spatial inhomogeneity, due to either 
engineered and disorder potentials, has been studied for 
both excitons [1| and polaritons 

By looking for steady state solutions to the cGPE, we 
find that a density-dependent gain rate combined with 
spatial inhomogeneity leads to steady-state currents, con- 
necting regions of net gain with those of net loss. These 
supercurrents in turn affect the density profile (as is al- 
ready well known in the case of solutions with vortic- 
ity), and so pumping and decay can significantly alter 
the density profile of a trapped condensate. The effects of 
steady-state current flows in the absence of pumping were 
considered in Refs. [25|, |26|. By studying the stability of 
these steady state solutions, one finds that with homoge- 
neous pumping these solutions become unstable to break- 
ing of rotational symmetry. Stability can be restored by 
reducing the size of the pumping spot to be comparable 
to the self-consistent size of the condensate cloud (set by 
the balance of pumping and decay). By increasing the 
pump spot size (or by decreasing the pump strength), 
the rotationally symmetric solutions again become un- 
stable, and are replaced by solutions with vortex lattices. 
The observation of vortices driven by the combination 
of particle flux and spatial inhomogeneity has been seen 
experimentally [5|; our results indicate that such vortex 
solutions can arise even with symmetric traps. Our find- 
ings show the existence of new phenomena in the already 
rich world of complex Ginzburg-Landau equations [22| 
that play an enormous role in our understanding of non- 
equilibrium physics and pattern formation [l^ • 

Our cGPE can be derived as the gradient expansion 
of the saddle point equation of a non-equilibrium path- 
integral theory of polariton condensation [19, 20]. How- 
ever, to provide insight into its meaning, we instead de- 
scribe here the physical origin of the terms it contains. 
The form of the cGPE depends on whether one consid- 
ers coherent or incoherent pumping. Coherent pump- 
ing, injecting particles directly into the condensate at an 
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energy cjq, is described by a source term 9tV^ = Fe*^^* 
[28^. We instead consider non-resonant pumping, and 
thus we introduce stimulated scattering into the conden- 
sate, 9t '01 gain = 70^ • A similar term ^t^|^\loss = —i^'^ 
describes particle decay, i.e. loss, and so we introduce 
7eff = 7 — With such gain and loss, the dynamics 
is unstable and trivial; if gain exceeds loss, the conden- 
sate grows indefinitely, if loss exceeds gain, the conden- 
sate vanishes. In practice, for non-resonantly pumped 
solid-state systems, the gain is saturable — it tries to 
bring the condensate density into chemical equilibrium 
with some external particle density. The simplest model 
of such a process is a density-dependent rate of gain, 
StPlgain = (7 — rp)p, which trics to establish equilibrium 
at p = 7/r. A closely related model of saturation, con- 
sidering a reservoir of non-condensed particles was stud- 
ied in [12]; the steady state behaviours of both models 
are very similar. We combine these terms and write the 
complex GPE in the following form 
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where V{r) is an external trapping potential, and U is the 
strength of the ^-function interaction (pseudo) potential. 

We will we look for steady state solutions and intro- 
duce the chemical potential, /i, in the usual way, via 
i?idtil^(t) = In this equation /i is a free parameter 

to be determined from the balance of gain and loss; nei- 
ther the chemical potential nor total number of particles 
is externally imposed. We will illustrate how the inter- 
action of spatial inhomogeneity with pumping and decay 
modifies the density profile by studying how the profile 
depends on pumping strength in a number of cases. 

We consider the classic example of a parabolic trap- 
ping potential in two dimensions. For this problem, two 
dimensionless parameters control the behaviour. We can 
write the potential as V{r) = (fia;/2)(r^/P), where co is 
the oscillator frequency and / = y/h/muj is the oscillator 
length. Expressing lengths in units of /, energies in units 
of huj, and rescaling y^huj/2Uilj, yields: 
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For the rest of this letter, we shall write ji = 2fi/huj^ 
and introduce the two dimensionless parameters which 
control the density profile: a = 27eff/fia;, and a = T/U. 

Before discussing the solutions, we give illustrative 
values Of, a, relevant to the polariton experiments of 
Refs. [3I, HJ]. The maximum pumping strengths con- 
sidered are around ten times the threshold pumping 
strength; this threshold occurs when pumping matches 
decay rate j = hz. The decay rate, found from the 
linewidth at low power, is :^ O.lSmeV, and so the pump 
rate may be up to 7eff < 1.2meV. To find a, one needs 



also the characteristic trap scale. In Ref. [2J, the dis- 
order traps are estimated to have a depth £^0 — 0.5meV 
and size a :^ 3/im, which with a polariton mass of m :^ 

10~^mo yields a trap frequency hu = \J Eoh^ /mo? 
0.2meV; hence 0<a^lO. Fis harder to estimate with- 
out a specific microscopic model; an order-of-magnitude 
estimate may be found from the observed blue shift (shift 
of chemical potential) vs pumping power. As discussed 
below, for weak pumping one has ji (ficc;/2)(3a/2cr), 
and so cr 37eff/2/i. In Ref. [3], a pump power at 
twice threshold, i.e. 7eff ^ 0.13meV yields a blue shift 
ji ~ 0.5meV, giving a 0.3, however this estimate in- 
volves considerable uncertainty. 

We first discuss the rotationally symmetric steady 
states of Eq. ([2|), using fixed point iterations combined 
with the secant algorithm for determining ji for a variety 
of parameters. We compare the densities of the ground 
state with the analytical Thomas-Fermi (TF) profiles 
found by neglecting density gradients and assuming that 
super currents do not affect the density distributions. Fig- 
ure [T] shows the density profiles for different values of a. 
As a increases, two effects are clear: firstly the increased 
pumping rate evidently leads to an increased total den- 
sity of the condensate; secondly increased pumping leads 
to a greater fiux, and so for a = 4.4, the density profile 
is not the TF profile, but is suppressed in the middle, 
where the super current is highest. The increase in total 
density can be described from the balance of net gain and 
loss; by multiplying Eq. ([2j) by V^* and integrating over 
all space, the imaginary part of this equation is: 
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When pumping is not too strong, substituting the 
Thomas- Fermi solution = {(1 — r'^) for r < into 
this equation yields /i = /i* = 3a/2cr. 

The suppression of density due to supercurrent means 
that with increasing pumping, the density profile be- 
comes increasingly sharp as supercurrent fiows become 
important. Such results have been strongly hinted at in 
several microcavity polariton experiments where sharp- 
ening of the peaks of the density profile with increasing 
density is seen with both disorder traps 2^, 2^ and en- 
gineered stress traps [4]. 

Let us now discuss how the changes to the density 
profile seen above can be understood physically. The 
Madelung transformation, = y^e*"^, represents Eq. ([2]) 
as a continuity equation and Bernoulli's equation: 
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Regions of high density imply loss, and regions of low 
density gain, which lead to supercurrents V(/), between 
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FIG. 1: The densities, p(r), of the steady states of Eq. (|2]) for 
a = 0.3 and a = 1.5,4.4 (black sohd hnes) as compared to 
the TF solutions p — fi" — r^^r < y/JI* (dashed lines). The 
inset shows the symmetric solutions for a — 4.4 with finite 
pump spot size as labelled; = 2,4 are stable, while R = 8 
is unstable to breaking rotational symmetry. 

these regions. If these supercurrents are large, they af- 
fect the Bernoulli equation, leading to a density depletion 
where current is largest. This is clearly seen in Fig. [TJ 
there is net gain at large radii, and net loss at small 
radii, and a dip in the density profile in between these 
indicates a region of maximum supercurrent. The radial 
phase gradient associated with the current could be seen 
experimentally using interferograms as in e.g. [1, [H, [l^; 
the maximum phase difference across the cloud scales as 
(X cr/i^ and 30, for relevant parameters. 

To study stability, the time evolution of Eq. ([T]) is fol- 
lowed, using the rotationally symmetric steady- state so- 
lutions as initial conditions, and including a small per- 
turbation. With an infinite homogeneous pump, as in 
Eq. ([1]), the solution is always unstable to angular per- 
turbations. This instability can also be seen by consid- 
ering pumping and decay as perturbative corrections to 
the hydrodynamic modes of a trapped 2D condensate; 
one finds the leading order correction to the mode ener- 
gies introduces growth/decay rates which always produce 
growth for modes with large enough angular momentum. 
Physically, this instability can be understood by look- 
ing at the region just outside the condensate cloud. In 
this region the steady state gain is zero since it is pro- 
portional to density, however linear stability analysis for 
p ^ p -\- Sp depends on dp [{a — (Jp)p] oc a — 2(jp, which 
is positive outside the condensate cloud, so any small 
perturbation will grow. High angular momentum hydro- 
dynamic modes of the condensate are unstable because 
they transfer density to the edge of the condensate. 

This mechanism of instability is supported by observ- 
ing that the instability is not present with a finite spot 
size; this ensures that outside the condensate cloud there 
is no gain, and so no growth. For simplicity, we treat this 
radial cutoff by replacing a by a{r) = a&{R — r), where 
9 is the unit step function and R the cutoff radius. A 



finite spot, of size comparable to the observed cloud is in 
fact used in current experiments [1, 0, • For small R^ 
this stabilises the radially symmetric modes. However, 
when R exceeds the Thomas- Fermi condensate radius, 
y/p y/3a/2a, the instability reappears. The subse- 
quent time dynamics, leading to a new steady state is 
shown in Fig. [2l The final state is no longer stationary, 
but instead rotates according to: ifidtil^ = {p — 2VtLz)il^^ 
where Lz = i{xdy — ydx)- 




FIG. 2: (Color online) Time evolution from the rotationally 
symmetric steady state solution of Eq. (|2]) for a = 0.3, a = 
4.4 when the radius of the finite pumping spot is = 5 (as 
marked by the red line). Times are in units of 2/uj, were uj is 
frequency of harmonic trap. 



As the initial problem is rotationally symmetric and 
non-rotating, the vortex solution spontaneously breaks 
rotation symmetry; either sign of vortex array is stable, 
but the rotationally symmetric solution is not stable. 
This behaviour is characteristically different from the 
equilibrium non-rotating trapped condensate in which 
vortex solutions are unstable [30], and vortices would 
spiral out of the condensate — with pumping and de- 
cay the dynamics shown in Fig. [2] shows that vortices 
spiral into the condensate. In addition, for a given ra- 
dius of pump spot, more than one vortex array may be 
stable, the number of vortices depending on the history 
of the spot size; this is indicated in Fig. [3l This is simi- 
lar to hysteresis effects in rotating Bose-Einstein conden- 
sates [311 , but the external rotation is absent in the model 
considered. The origin of the instability — growth of 
condensate density outside the Thomas-Fermi radius — 
suggests that other models of the cGPE with a reservoir 
would show the same behaviour [18]; for the instability 
to be removed, one requires a reservoir concentrated near 
the minimum of the trap. In the context of the polariton 
condensate, this means the instability might be cured if 
the non-condensed exciton reservoir was highly mobile 
— in the language of laser theory, this corresponds to 
damping of instabilities by carrier diffusion. 

This solution of the cGPE can be understood as vor- 
tices enlarging the cloud size to match the pump spot. 
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Adapting Eqs. ()4|5p for a rotating solution [2l| gives: 



V • [p(V0 -nxr)] = {ae{R - r) - ap) p, (6) 



(7) 



The rotating vortex lattice solution adopted can be un- 
derstood as follows; vortices lead to quantised rotation, 
and the density of vortices, o:^ Vt/n ensures that 
Vt X V mimicking solid body rotation. Neglecting 
the vortex core, the continuity equation, Eq. (|6|), thus 
requires p a/ a. This implies that inside the vortex 
lattice there is no net radial current in contrast to the so- 
lutions with smaller R; for the solutions with a single vor- 
tex (when R < ^/3a/2a)^ the combination of radial and 
rotational currents means such vortices are in fact "spiral 
vortices". For this constant density solution to be valid 
(except near each vortex core) Eq. ([7j) requires 1] :^ 1 and 
p = p a/ (J. This solution persists till the edge of the 
vortex lattice, beyond which = Ny/r, where A^v is the 
total number of vortices. The total number of vortices 
is then set by requiring the edge of the cloud to occur 
around r = leading to n^jnR^ :^ R^ for large R. 

When the vortex core is not negligible, the extra gain in 
the vortex core and quantum pressure corrections imply 
p > p > a/ a. This is shown quantitatively in Fig. [3l 




FIG. 3: (Color online) Chemical potential /i as a function of 
pumping spot size R. Numbers of vortices for the stable vor- 
tex lattice configuration are marked. Solid lines join branches 
of stable solutions. Dotted line corresponds to unstable radi- 
ally symmetrical solutions without vortices. Horizontal solid 
line marks jl — a /a. Inset: vortex lattice for = 10. 

For yet larger R^ e.g. = 20, no such simple rotating 
vortex lattice is found — for such parameters there is a 
residual vortex lattice in the center of the cloud, but the 
behaviour at the edge becomes irregular. 

In conclusion, we have shown that steady-state cur- 
rents connecting regions of net gain and loss can lead to 
significant modifications of the density profile of a quan- 
tum condensate, leading even to instability of the rota- 
tionally symmetric state and the spontaneous creation 



of a vortex array. Vortices can be clearly observed as 
in Ref. [5[ by pairs of forks in the interferogram of the 
emitted light. To observe the spontaneous vortex array, 
one would require that disorder is weak compared to the 
harmonic trap, which may prevent its observation in the 
current generation of semiconductor microcavities, how- 
ever other than this hurdle, the numerical estimates place 
current polariton experiments in a regime in which such 
effects could occur. 
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